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TFCHNICAL MKMORANDUM 7*<142 


SINGLE NODE ORBIT ANALYSIS WITH RADIATION 
HEAT TRANSFER ONLY 

INTRODUCTION 


In many caaes, it is di-sirable to simulato an orbital analysis of a 
thermal resixjnsive system without applying I.OHAHP or SINDA. These are 
library programs capable of handling the most complex and sophisticated 
thermal system. However, in many cases, tlie depth which can lie pentirated 
by these programs are not warr;inted on the basis of required turnaround time 
and technical level, 

A need exists, tliereforc, to handle an orbital thermal problem where 
the fidelity of the analysis is commensurate with th<- level of decision making 
demanded by a phase A/D study. This includes problem modeling and problem 
solution to be accomplished within a few days. 

One approach to satisfying tliis need is to consider seh-cted types of 
thermal systems. The purpose of this report is to provide background tech- 
niques leading to orbital analysis of a single node having radiation heat transfer 
only. To this end, a computer program has been developed to handle this 
category of thermal systems. Specifically, these* arc; 

1, Ix)uvers systems 

2, Flat plates. 

These two basic categories are sulxiivided into five combinations involving 
thermal systems having shields, sunshades, and louvers. Each combimition 
is given further definition in ’Thermal Transients Resulting from a Time 
Varying Thermal Environment," contained herein. 

For illustration, tlie same inputs have l»en applied to each combination. 
Unless otherwise stated, these conditions are: 

1, Orbit period — 90 min 


2, Orbit day period — 60 min 
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Solar conHtant — -130 IUu/l»-ft* 


4 , Albedo and IH — 30 lUu4i-ft* (orbit aver:me) 

5. r>ptical proiK*rtio8 — 

a. Kquipim-nt eniittanci* — 0,80 

b. Shield einittunce — 0, 80 

c. Sun shield eni ittanee — 0, ho 

d. Sun shield ab8oi*))tion — 0. 10 
G, Kquipment capacitance — 3 Htu/*F 
7. Initial nock* temperature — 70»F. 

These inputs have been selected for example punioses only (no s)nniflcance is 
to be attached to them). All are variables within the proRram and can lx* 
changed to suit a specific problem. 


STEADY-STATE TEMPERATURE CHARACTERISTICS 
RESULTING FROM A NON-TIME-VARYING 
THERMAL ENVIRONMENT 


In considering steady-state results, it is im|)ortant to differentiate 
between a time-varying Uiermal environment and a fixed thermal environment. 

In general, the results are not the same. In the general case, the thermal 
capacitance of the system prevents the temiwrature from increasing or decreasing 
to the level prc'dicted by a fixed thermal environment. Sometimes, the parame- 
tc'rs of the sysU'm can be such that very grc*at differences exist. 

Consider a flat plate c'xposed to a radiation flux, C'l, as illustrated in 
Figure 1. The plate has pro[M?rtie8 as indicated with insulation on the antiflux 
side. An internal heat source, q/A, i*epresents the presence of electronic 
equipment that may be attached to the plate. In the steady-state, all of tlie 
energy in|)Uts to the plate must lx- equal to the energy radiated; 


I 


I 




FJKurt* 1. Flat plate showing enc*rny i.*omponentH, (The entire 
plate is considered as a single node.) 


oG + ^ = a e T ^ 
A 0 0 


(0 


Solvint; for tempo r.iture, T^, 


T « 
o 




( 2 ) 


The sU’.idy-stato temperature as defined in equation (2) is indi>pendent of plate 
mass. The plate is at a uniform temperature with no pradic-nts. Actually, 
equation (2) is a stronp tool for evaluatinp the combined effect of optical pro- 
perties, internal heat peneration, and incident radiation. However, unless the 
system has low mass per unit aix'a, additiomil tuialyses are warranted. 

Note that the plate temperature and surfcice emittance has been sub- 
scribed with the letter (o). This subscript is meant to imply the surface on 
which electronic equipment is mounted. This clarification is helpful when 
shields are considered. 
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SHIELDS 


The next level of sofihistication to lx* considered in this tyjx problem is 
the introduction of a shield. Sometimes, shields are referred to as sunshades. 
This is especially true if the shield is used to inU‘rrupt solar radiation from 
impendinK upon the plate. WiUi the use of multishields, the outermost shield is 
referred to as the sunshade. 

The Beometry of a shield is illustrated in Figure 2. The shield is placed 
in fro ‘ of the plate. The outer area bi‘tween tlx shield and plaU* is di'siBned to 
be as small as possible. For purposes of analysis, this art*a is assumed to lx* a 
thermal barrier. 



Fipure 2. Shield configuration with surface 

emissivitv of <r . 

n 

In application to thermal control systems, the shield concept has the 
ability to limit the temjx'rature extremes predicted by a single plate. The 
thermal cycle is ’’smoothed out.” 

The steady plate temix*raturc, T , can be obtained by writing an enersy 

o 

balance Ix'tween the shield and plate. Under steady-stati- conditions, all of the 
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equipment energy mu«t be radiated to the ahield. ThuH, the 8teady-Htati' tem- 
perature of the shield from equation (2) Ix'comes 


'^shield 



( 3 ) 


The energy balanct* betwi*en the shield and plate is 


a AIT * - T<J 
o 



( 4 ) 


Substituting equation (3) into equation (4) yields 


a 

A 


U 



1 




It., a T * -ofi) . , . 
rs o 


( 5 ) 


Equation (o) can be easily solved for T , l)ut is much less complicated in this 
form, ° 


It is noted tliat equation (5) applies for a single shield having an 
cmissivity value of e., on botli outside and inside. If the inside value is e, 

equation (5) becomes 


a 

A 


rs 


lCj^,aTo«-ar.l 



1 


(G) 


This equation is important Ix^cause it illustrati’s how the heat dissipation 
capability can be reduced by making c a small value. In ceicain applications 
it may be desirable to reduce tlie heat dissapation rate. Also, a small power 
input may be desirable to maintain a constant temperature. The rate loss can 
be easily reduced by a factor of ten by adjusting the value of c. 



Analynis of Hhifldr and 8unHhack.-8 arc directly a|>|*Ucable tu louver Hyn- 
tcniH. A louver Hyotem iH one where the Hurfaoe area Ix'int; Hhuded is variable. 
Thus, tlte rate of heat dissipation can lx* varied over a wide ranne. (ieni'rully, 
thi're aix> twt) louver conf);oiratiuns of interest. These are Ivnuwn as "wiU> 
sunshade" :uid "witliout sunshade," A louver sy8t(*ni uitli sunshade is known as 
an "inU'rnal louver system," 

For the intL'mil confitnii'ntion, one shield would represent a "full^ixm" 
louver situation, and two shields would represent the "full*closed" situation, 
Equation (0) represents Ute full- 0 |K‘n situation. This equation is plotU d in 
Figure 3. Fmissivities have bet>n selected to give ni‘ar maximum heal dissipa- 
tion capability. Equat ion (2) represents the full o|ien situation for i^e "without 
sunshade" configuration. 

The same U'chniques as already presented can be a|>|)lied to two shields 
as illustrated in Figure 4. This con^ig.iration is a|)plicable to an inU>rn:il louver 
system with the louvers closed, I’nder sU-ady-state conditions the equipment 
ent'fgy which can be radiated for the configuration illustraU-d in Figure 4 is 


A “ S' 


T * 

+ c . c . 

o 

2 3 

C, 

c„ + 1 

1 

2 


- n ( '. 


Cj » l2-«) 


c+3c -2c c 
o o 


c (if - t) 
o' ' 


["f-1 V 


(V 


Equation (7) is plotted in Figure 3 for typical values. The |x)sition of 
the two lines can lx* changed by varying the emisslvities. However, the position 
shown is for a value of 0,8. The opened louver position curve represents thi> 
near ouU*r capability of louvers to control temperature. 


G 
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FiRurc 3. Steady-staU pcrformanco characteristics of a louver system 
in a nonvaryins thermal environirent. 







r» ,urv •!. ConfiKurutiun of a flat plaU- with 
Hhk'ld and Hunshndo, 


It |8 important to recognize tiiat Figure 3 applies to steady-state condi- 
tions on the day sick* of Uie orbit. An itientieal curve can lx* develoix'd for the 
night cycle (Fig. 5). However, as will be illustrated later, Figure 5 is a j)oor 
compromise for results of orbit:d lUialysis. As a consequence, day conditions 
are sometimes substituted for orbital analysis results. It will lx* demonstrated 
later tliat orbibil :uialysis of louver systi*m ivsults in gix*aU*r heat dissipation 
thr^' indicaU'd in Figure 3. 


TEMPERATURE TRANSIENTS RESULTING FROM A 
NON-TIME-VARYING THERMAL ENVIRONMENT 

Knowledge of tlie transients occurring \\ithin a thermal system adds to 
the understanding of the problem and inspires solution techniques. The approach 
to establishing the transient can be state<i in the form of an instantaneous heat 
balance: 


M C 


p dr 



(«) 
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whore is the not onorKy rato fjaln (or loss) of the system. In general, 

JjQ is represented oy a radiation term, thermal flux, and equipment power. 
If the equipment power :uid thermal flux are zero, equation (H) reduces U> 


dT 

M C T- - 
p dr 


This expression can Ik‘ easily integratiHi to obtain a tempi'rature time expression. 
For those systems involving shields, heat flux, equipment lower, and u time- 
varying environment, a computer analyzer progiam is the only practical way to 
approach the problem. However, there is a special case tiiat can ixi h:indied 
mathomatically. This is u simple flat plate under radiation heating, not subject 
to a varying thermal environment. Kquation (H) becomes 


dT a c \ T* o G A 

dT M C " M C 
P P 


( 10 ) 


Solution to this equation can lx* found in the literature, but will be repeated here 
for pur|>OHe8 of completeness: 



T^ = initial plate temperature (*H) . 

Equation (11) gives temperature as the independent variable. However, tliis 

jl4 

equation can bo generalized upon by plotting t ~ as an independent variable 

(Fig. 6). Note that in a mathematical sense, steady-state temperature is 
obtained only after an infinite time. The value of the steady-state temperature 
can be determined by 
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T (steady-state) - 


The argument of this parametric plot is b'a, reciprocal of the steady-state 
temperature (Rankine). The independent variable allows representation of all 
combinations. 


T — » T 

a 


However, Figure G is valid only for an initial plate lem|K‘rature of 70*F. Thus, 
a complete generidization of equation (ll) cannot lie made. 

If it is desirable to determine just the time to reach the steady-state 
temperature, it becomes necessary to talk in terms of "time to reach a per- 
centage of steady-state* tem|K*rature," This is necessary since in a mathematical 
sense an infinite time is rwjuired. To solve for the time to achieve 0.9H percent 
of tile steady-state teni|x*rature, let 


T = O.QS^ » '***^'*[ 7 ^] 


In most cases tlie tangent functions can be eliminated since the equation is 
dominated by the logaritlim function. This fact allows for quick calculations. 

The following is a sami?le case to use in conjunction with Figure <>; 

How long will it take a solar array to reach ISO'F under the following 
conditions: 


T^ = 70*F (initial temperature) 

1 / 

e = 0.8 , C = 0.2Btu/lb-»F , ^ = 2 lb ''ft^ ? 

D A 
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KntcT Finure (5 at 1M()*F for a (a/b) valuo of 0.0012: 


• 

P 


Making tht* nccoK8:iry MubKtituUonrf and Holvinu for r KivcK; 
T - 0.1770 - O.OSSS h (r..:i2 min) . 

U • K 


Tho v;iluo of IhiH t>pc c:ilc’ulation i8 oHtablishcKl if hufficicnt time exiHtK 
for steady-state e>>nditions. If time is not available, then a more sophisticated 
solution is \^:irranted. 

Note: how the emiUance has a dual effect. It affects the v;due of a^ and 
T a*; however, its effects on t a^ is much more sensitive. A j;reater emissivity 
will decrease the heiiting time. Also, the mass |x*r unit area is impoiiant. 

Kquation (11) h:is bt?en demonstrated as considerinj? a heating problem. 
Cooling problems can lx* solved by the same equation. For example, if T > 
b a, coolin^r will occur. Also, T must be less than b/a. 

As a matter of clarification of tin* units involved in i*quation (0), it is 
noted that 



must have units of time. At first this is not obvious. Hy making the substitu- 
tions into this set of parameters, the time unit of t can lx* verified: 
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- Ih-CK)*) - h . 

TEMPERATURE TRANSIENTS RESULTING FROM A 
TIME-VARYING THERMAL ENVIRONMENT 


The application of equations prescntcti prcviouKly for a fixed environment 
have limite<l utility. They apply only for a uniqu(‘ Het of circuniKtaneeH which 
the cnj^ineer seldom has the luxury to afford. 

At best, application of the previous equations represents no more tlian 
’•first impressions” and, in Kor.eral, caniKit nuH*t tlie feasibility definition 
required for a Phase A Level Study. As in exiuiiple of what is implied by this 
disposition, consider tlie tc'inperature of a flat plate in Kartli orbit having the 
following characteristics; 

f c 0. r> 

o 


o * b. 5 


a 

A 


W ft‘ 


From equation (1), the steady-state tempi*rature on the Sun side is 
2iil*F. On the ninht cycle, the steady-state temperature is ♦28* F. The actual 
temperature profile resulting from computerizi*d luialysis is illustrated in 
Figure 7. In the steady-state, the plate temperature will vary from a maximum 
of 212*F on the day cycle to a low of 19(»*F on the night cycle. The capacitance 
of the systems combined with the cycle cnvironnu'nt prevents Uie system from 
reaching peaks jiredictcd by equation (1). Note that steady-stati* conditions are 
achieved after 900 min (10 orbits). 

The re;il value of treating a proldem alxive the sophistication of equation 
(1) is the inlierent expansion of knowledge ai»ut the Ix'havlor of thermal sys- 
tems. The data of Figure 7 were developed from the computer program descrilx;d 
in tlie appendices. 
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0 200 400 600 800 1000 

TIME IN ORBIT (mm) 

Figure 7, Temix-Taturi* res|)onse of a flat plate*. 


THERMAL CONTROL ALGORITHM FOR LOUVER SYSTEMS 

One consideration that uxpamls knowledjfe of a louver system is the rela- 
tionship between temperature and the defleclion an^le of the louvers. This lela- 
tionship is called tlie thermal control iili'oritlim. This control law determines 
the actu:il equipment temperatures. The steady-state full-open closed charac- 
teristics given in Figures 3 and 5 only establish the other temix-rature boundaries 
that are tliormally possible for a given equipment ix)wer. All of these possible 
temperatures are not desirable. Figures 3 and 3 can lx? misleading since they 
do not illustrate actual temperatures achieved. Thermal capacitimce combined 
writh the thermal control algoritlim allows any of the temperatures intiicated 
(within small fluctuations) to be maintained during the entire orbit. 
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For pur|M)HL‘S of illuHtratinu how U!ni|x>rature can Ijc controlled, two 
tt-.ermal control alKoritlims will !*.• introduced. Tht-Hc are Klven in Figure K and 
labeled ( 1) :in (U). Ik>th have a slope of .'1 degrecH '*F, However, for(l), the 
louvers do n- open fully (‘JO degrees) until W)*F is reached. For (2), the 
louvers are laiiy open at 00* F. Thus, it is obvious that control law (2) will 
maintain tlie equipment at a lower U*m|K'rature. Figure 0 shows the steady- 
state orbital tem|M-rature for algorithms ( 1). The louver system accomplishes 
a near constant H0*F over Uie entire |»ower range of interest. Suporimjiosed 
upon this figure is the day side analysis of Figure U. Notice how trivial the day 
side analysis is compared to a simulated orbital analysis. 



NODE TEMPERATURE (°F) 


Figure H. Illustration of the two example Uiermal control algorithms. 

(Specific temperatures encountertxi by a louver system are 
determined primarily l>y the algoritlim employed). 

Figure 10 is equivalent information using algoritlim ( 2). Here, two cases 
arc shown with and without sunshade. As exjH'cted, this algorithm maintains a 
lower equipment temperature. The capability of a louver system c:ui lx- doubled 
by not employing a sunshade. I'se of a sunshade will depi'nd upfin the equipment 
power level and the range through which it must lx* modulated. Figure 11 illus- 
trates the temperature ;uid louver angle variation during the orbit. The louvers 
vary through an S degree :inglc while the trough is approximately a 1.5 degree 
change. 
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Figure 9. Steady-state orbit tt^mperatures resulting from :i|)plicutinn 
of algoritlim (1). (The p«.>rformance sensitivity to a sunshade is 
dramatic. In general, heat rejection capability- can lx* 
doubled by removal of the Sun shield.) 

It is noted that application of a shield increases the thermal radiation 
resistance. A system \\1thout sunshade allows for greatijr heat dissipation 
since there is less thermal resistanc-e to radiation. 

For purposes of providing the capability for so|)histicated analysis and, 
at the same time, one simple enough for quick turn-around, a computer program 
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Figure 10. Steady-state orbit temperature resulting from algorithm (2). 






TEMPERATURE 


SBiuoao-NOiioanjaa bSAnoT 



19 


Figure 11. Equipment tcmperuture response during on-orbit for equipment poue*r of 
16 W ft’. Louver :ingle movements arc also given. Only slight temperature 
variaMonr are encountered. The data resulted from 
application of algorithm (2). 




hns bvi'n dovelopod (or HtudylnK thv tvmpvraturu tranMirntu of a ainKk* nrxk* in a 
time varyinit thornml environment which conNidora radiation heat transfer or^y. 

A listinR of thin program and an explaration of die limitutionH and constraints 
are Riven in Appendix A. Appendix n develops the rationale and assumptions used 
in this program. Tho computer program is flexiiile in that five types of thermal 
systems can be facilitated: 

1. liouver system with sunsh: Je. 

2. Ix>uver system without sunshade. 

3. Flat plate without louvers or sunshade. 

■1. Flat plate with one fixed shield. 

5. Flat pl:ttc with two fixed shields. 

Comments stated within the program indicate how each tiicimal system may 
b<? executed. 
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APPENDIX A. LIMITATIONS, CONSTRAINTS. 
AND OPERATION OF THE COMPUTER 
PROGRAM 


In applyiiiK th<- cunjpuUT thr o|KTutor Hhould Ik* ;iw:irt* of lh<* 

followInK UmltatlonH and inxtrut'tionN. KxctcIho Judt'c*nu*nt in ili*v»*I«)pnu*nt of 
thi* inputs. This i-autlon in c-on8i8tc*nt uith tin* Kimplicity rcprcKontcd in i*M*cut- 
inu tht* program for ((uick model development atxi turnaround. 

1. CommentM arc provided U> eNecute om of the five conficurntions 
nvnilahle. The instruetionH cauHe Uu* program to default to the pro|K*r confiKU- 
ration. Note* the H|M*eific8 of the followinc KtatementK; 

a. Stat«.-ment 71 — Orbit time (min). 

b. Statement 75 — Sun time ( min). 

c. Statement 7(5 — View factor, Karth tt) nfxle. 

d. StaU'ment 91 — Solar ant* iillH*do flux (Htu/h-fP). 

e. StaU-ment s7 — Kartli IH ( Mtu h-ft‘). 

f. Statement lUx — ^ives tem|M*rature at a 90 denr<*e louver angle. 

g. Statement 109 — giveH temnt*rature at a 0 degree louver angle. 

h. StaU*ment 110 — Algorithm, equation ol louver ;mgle versus 

tein|K*rature. 

i. Statement 95 - ixode view factor U» space on day cycle. 

j. StaU*ment 97 — Node view factor to space on night cycle. 

2. The program, as presented, d(H*s not accommodate a Hat plate with 
a]lx*do. As written, the flat plate is insulated mi one side. For a flat plate 

w ithout insulation, replace equation in statement 19 with QT » IF(93) + F(97) |* 
EN^SIG* I) -G*AC. 

3. Start point is at the terminator toward the day side, 

4. The time into and out of the terminattir is considered U> lx- instan- 
tam'ous. Statem<*nt H9 determines the day side flux and statement h 5 dc*termim s 
the night side flux (IH). 
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APPENDIX B. DEVELOPMENT OF RATIONALE FOR 
ORBITAL ANALYSIS OF A SINGLE NODE EXHIBITING 
HEAT TRANSFER BY RADIATION ONLY 

Consider the flat surface shown in Figure H-1 which him external heat 
{generation distributed equally over the surface. In the Kener:d case the surface 
is shaded by a louver and sunshade. The louver rotated through an^le 0 which 
varies the area shaded by tlie louver. At 0 » IM) de{;ree, the flat surface is 
shielded only by the sunshade. At 0 » 0 denree the flat surface is shaded by 
two shields. 

Mathematically, the temperature of tlie plate with ont‘ effective shield is 

T., and the temperature of the plate witli two effective shields is T_. The bulk 
N C 

temperature, T, of the plate becomes 


Me ♦ \ 


Me * M^ 


(B-1) 


where M_ and M are the respective shielded mass of the flat plate. The bulk 
C N 

temperature, T, is defined as the node temperatui'e and is the combined effective 
temperature of the shielded plate. 

It cim Ijo easily shown that and the projected area, A^ is 


A^ = A cos 0 


(n-2) 





cos 0 


(B-3) 


Me + = M 


(n-1) 
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SOLAR FLUX 



FiRuic H-1. General configuration of an 
internal louver system. 


and 

= A (1 - cos n) 




+ M 


N 


1 - cos 0 


where A is the total area. 


The equation describin^j becomes 


dT / V 


ft) ■ ft) 

\ / r*r\irw rfc\ \ ^ ' 


(EQllP) (RADIATED) 


C 

N p dr 


a Aj^,G(t) + 




'(EQUIP) 


(RADIATED) 


(R-5) 

(R-G) 

(B-7) 

(B-H) 
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Those t>»‘o equations can be represented as a single equation involving only the 
node temperature, T; 


A p 




(HADIATKD) 


(H-9) 


The term (q/A) equipment is simply ;ui in|)Ut equivalent to the ecfuipment |x>wer 
applied to the node, Kquations ((’») and (7) of ti!Xt can lx* combined to obtain tlie 
equipment power capable of Ix'inn radiated. It is noted that the solar flux falling 
on surface A is a function of time dependinn upon orl)it par:uiieters and orienta- 
tion. 


Kquation (B-9) is the problem representation included in Appentiix A. 
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